The steady state quantum statistics of a non-Markovian atom laser 



J.J. Hope 1 ' 2 '*, G.M. Moy 2 , M.J. Collett 1 and CM. Savage 2 
Department of Physics, University of Auckland, New Zealand 
2 Department of Physics and Theoretical Physics, Australian National University, ACT 0200, Australia. 

* email: jjh@phy.auckland.ac.nz 
(February 1, 2008) 



0\ 



3 

(N 



> 

m 

(N 
O 

o 

OS 

1 



c 

o 
o 



X 



We present a fully quantum mechanical treatment of a 
single-mode atomic cavity with a pumping mechanism and 
an output coupling to a continuum of external modes. This 
system is a schematic description of an atom laser. In the 
dilute limit where atom-atom interactions are negligible, we 
have been able to solve this model without making the Born 
and Markov approximations. When coupling into free space, 
it is shown that for reasonable parameters there is a bound 
state which does not disperse, which means that there is no 
steady state. This bound state does not exist when gravity is 
included, and in that case the system reaches a steady state. 
We develop equations of motion for the two-time correlation 
in the presence of pumping and gravity in the output modes. 
We then calculate the steady-state output energy flux from 
the laser. 
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I. INTRODUCTION 

Since the experimental realisation of a Bosc Einstein 
condensate (BEC) in a weakly interacting gas there 
has been a lot of interest in using similar techniques to 
produce a superior source for atom optics experiments, 
which are limited by the linewidth of thermal atomic 
sources. It was seen that to produce an atom laser, it 
was simply necessary to generate a BEC and then co- 
herently couple it to the outside world ^ E3 . This was 
first achieved with sodium atoms at MIT by coupling a 
BEC from a magnetically trapped state to an untrapped 
state using an rf pulse |lq| , and has since been repeated 
with long rf pulses [jl6| and Raman transitions jlTj . Al- 
though the resultant matter waves have been the most 
monoenergetic source of atoms that have yet been pro- 
duced, a gain-narrowed atom laser produced with a con- 
tinuous pumping mechanism will have a spectral density 
which is orders of magnitude larger. In this paper we 
will describe a fully quantum mechanical model for such 
an atom laser which does not make the (invalid) Born 
or Markov approximations, and also does not make the 
mean field approximation. 

While a condensate is held in a trap, it is in an eigen- 
state of the system, and is therefore completely monoen- 
ergetic. If the trap is suddenly turned off or the conden- 
sate is quickly coupled into free space, then the result- 
ing wavepacket will have a spread in energies due to the 



momentum spread of the trapped wavefunction. Such 
a coupling tends to conserve momentum. Coupling the 
atoms out slowly will tend to preserve the energy of the 
intial state, so a monoenergetic output will be achieved, 
but this is at the cost of reducing the atomic flux M [l8| . 
A continuously pumped laser can have the best of both 
worlds. Due to a competition between the pumping and 
the damping, it can produce an increasingly narrow en- 
ergy spectrum in the output as the pumping and the 
output flux are increased. It is this feature that we wish 
to discover in our atom laser model. Several theoretical 
models of a continuously pumped atom laser have been 
produced |fl|-|l4||, but they have all made the Born and 
Markov approximations for the output coupling, which 
have since been shown to be invalid unless the output 
coupling rate is extremely small | IS - 20 1 . In this paper 



we produce and solve a model for a continuously pumped 
atom laser which does not make these approximations. 

The early atom laser models are largely distinguished 
by their choice of cooling method, which was either some 
form of optical cooling ||f7]||], or evaporative cooling 
Hf|. Evaporative cooling appears to be less appropri- 
ate for a continuous process, but it is the only method 
which has experimentally reached the quantum degener- 
ate regime in BEC experiments. In all of these schemes, 
the model for the damping of the cavity was the same 
as that used in the master equation for the optical laser. 
The resulting equations were therefore very similar to 
optical laser equations. This means that they could be 
solved using similar methods, and were shown to produce 
analogous behaviour. We show here that a correct de- 
scription of the output coupling leads to an irreducibly 
non-Markovian damping, and can lead to different be- 
haviour. 

In an attempt to produce a more realistic description 
of the output coupling, we modelled a cavity from which 
atoms were coupled to free space via Raman transitions 
JLlj . This allowed high intensities, spatial control and the 
possibility of giving the atoms a momentum kick from 
the lasers, which gave the outcoupled beam a direction. 
This was recently achieved by Hagley et al. Jl7| . We also 
suggested placing the beam in an atomic waveguide such 
as a hollow optical fibre, which is a possible method of 
achieving good spatial properties, and an effectively one 
dimensional output. Although the rate equations and 
the practicalities of the Raman coupling scheme appeared 
favourable, it was necessary to produce a more complete 
theory in order to describe its full effects. 
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A quantum mechanical theory for the output coupling 
from an atomic cavity connected to an external field was 
then developed . This theory described the dynamics 
of a BEC which is coupled to the outside world, but is 
not pumped by some continuous process. It was based 
on the optical input-output formalism developed by Gar- 
diner and Collett , but it was complicated by the dis- 
persive nature of the atomic field as compared to the 
linear energy spectrum for optical fields in free space. A 
general solution for the dynamics of the cavity and the 
external fields was presented in terms of Laplace trans- 
forms, and an analytical solution was found for a simple 
form of the coupling [|l9| . 

The quadratic nature of the dispersion relations in free 
space make the dynamics of atom-optical output coupling 
very similar to those of photon emission in materials with 
a photonic band gap. Near the edges of these bands, the 
dispersion relations are quadratic rather than linear, and 
this causes the model to behave in a non-Markovian way. 
The same qualitative behaviour was found in treatments 
of this system as we found in our purely damped case for 
the atom laser p2| , p3[ . 

A complete model of an atom laser must also include 
the effect of interatomic interactions and a pumping 
mechanism. The effects of atom-atom interactions are 
included implicitly in the output coupling models which 
describe the fields using the nonlinear Schrodinger equa- 
tion (NLSE), which includes the effect of s- wave scatter- 
ing (24|-|27]]. These models based on the NLSE do not con- 
sider the effect of the coherences between the lasing mode 
and the external modes, and therefore implictly made a 
Born approximation. In fact, all mean-field models are 
essentially semiclassical models, as they assume that the 
atoms in the lasing mode can be described by a spa- 
tial wavefunction with a predetermined set of quantum 
statistics. This makes it impossible to fully describe the 
behaviour of the pumping, and the dephasing of the laser 
mode must be added a posteriori. The full quantum me- 
chanical description allows us to calculate the quantum 
statistics using a microscopic model and physical param- 
eters, and therefore allows us to calculate the linewidth 
of the resulting output. 

Unfortunately, interatomic interactions are very diffi- 
cult to include in a full quantum mechanical model. We 
present a model of an atom laser which includes a pump- 
ing mechanism but assumes that atom- atom interactions 
are negligible. By making this approximation we have 
been able to model the output coupling without mak- 
ing the Born-Markov approximation, which is invalid in 
most physicially relevant parameter regimes. This is an 
accurate description of very dilute systems, and we show 
that it can be made self consistently with a laser oper- 
ating well above threshold, as the threshold can be less 
than 10 2 atoms in the laser mode |p8[ . This restriction 
is made for calculational purposes rather than as a delib- 
erate design criterion. The inclusion of the atom-atom 
interactions in full generality would require a multimodc 
description of the intracavity field. As discussed above, 



the mean-field methods which have been so successful in 
simplifying this procedure cannot be used without de- 
stroying the very information that we are trying to find. 

It may even be possible to use the single-mode ap- 
proximation for the cavity after the system has reached 
a steady state, if this mode is derived self consistently. 
At this time, the number of atoms in the cavity may be 
reasonably well defined, and the complicated dynamics 
of the pumping process may be approximated by a lin- 
earised master equation term, as may be constructed for 
an optical laser p9[ . 

Section || describes our model of a pumped and 
damped single-mode atomic cavity, and descibes the 
method by which we may determine the properties of 
the output field if we k now the dynamics of the mode in 
the cavity. Section III shows how the solution behaves in 
the absence of pumping. Section |y] derives the equations 
of motion for the system in the presence of pumping and 
discusses the physics of the various limits of the equa- 
tions. We see that we can write the equations of motion 
for the intracavity field in terms of the intracavity opera- 
tors only, but that this leads to non-Markovian equations. 
In Sec. |V|, we show how these equations can be solved to 
find the two-time correlation of the lasing mode. Section 
VI shows how the solution in the absence of gravity can 
be found largely analytically, and that it i s not self con- 
sistent as a steady state solution. Section VII describes 
the features of the energy spectrum of the outp ut of an 
atom laser in the presence of gravity. In Sec. 
discuss the possibilities for extending this model. 



II. THE MODEL 

We model the atom laser by separating it into three 
parts. The lasing mode is an atomic cavity with large 
energy spacing, and when it is operating in the quantum 
degenerate regime, it is effectively single-mode [|[^4). We 
assume that the cavity is single-mode, with annihilation 
(creation) operator and a Hamiltonian H s . The ex- 
ternal field has a different electronic state from the atoms, 
so the atoms are no longer necessarily affected by the 
trapping potential. We model the external modes with 
the field operators tp^(x) and the Hamiltonian H Q . The 
operators a and ?p(x) satisfy the normal boson commu- 
tation relations. The coupling between the lasing mode 
and the output modes will be described by the Hamilto- 
nian Hi. The pump reservoir is coupled to the cavity by 
an effectively irreversible process. At this stage, we will 
describe the pump by the Hamiltonian H pi which also 
couples the atoms from a pump reservoir into the system 
mode. 

The total Hamiltonian is then written 



VIII we 



H p + H s + Hi + H a 



(1) 



where 
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H s = hu) a'a, (2) 
Hj = ih I dx {n{x,t) ip(x)a' — K(x,t)* ip(x)^a) (3) 



and where 

[H ,H M ] = [i? ,a (t) ] = [H M ,^{x)] = 0, (4) 

as the pump does not directly couple to the external 
modes. This is described in Fig. 0. 



lasing mode 




FIG. 1. Schematic model of the atom laser 

We then enter the interaction picture, leaving only Hi. 
This gives us the interaction Hamiltonian: 

/oo 
dx {n(x,t) ip I (x,t)a\(t) - 
-oo 

K(x,t)* Mx,ty ai (t)) (5) 

where 

ipi(x,t) = e iHo{t-to)/n ^ Xih ) e -iH o [t-t )/n 



ai(t) 



_ i/H(H,+H p )(t-t ) 



a(t ) e 



-i/h{H s +H v )(t-t ) 



are the system and output field operators in the interac- 
tion picture. 

We therefore obtain 



H n 



-ih 



(m 4(*)-s f (*) <*/(*) 



(6) 



where = f dx K(x,t)tpi(x,t). 

Appendix H shows how we can go to the Heisenberg 
picture, and find the result 



il> H {x,t) = ipi(x,t) - / ds F(x,t - s)a H {s) (7) 

Jto 

where TpH(x,t) and an(t) are the Heisenberg operators 
corresponding to ip( x ) an d a respectively, and where 

F(x,t,s) - 

dy K*(y,s)G(x,t,y,s) (8) 



where G(x, t, y, s) is the Greens function propagator due 
to the output Hamiltonian, H Ql only. These functions can 
be written in closed form for several useful cases includ- 
ing: free space, free space with gravity, and a repulsive 
Gaussian potential. We may use Eq. [j] to calculate any 
observable in the output field that we desire, providing 
we know the complete history of the system, au{s). This 
is a simplified version of the input-output relations for an 
atomic cavity. 



A. Output field energy spectrum from (<r (t)a(t')) 

When the system is in a steady state, we are not in- 
terested in the external spectrum directly, as it is always 
growing. The quantity of interest is the output energy 
flux. We transform our interaction Hamiltonian into the 
basis of the energy eigenstates of the output modes: 

H = J dp h LJpclcp (9) 

Hi = ih / dp t) c p a* — R(p, t)* c], a) (10) 



where c p is the annihilation operator associated with 
the eigenstate of H D that has a position space wave- 
function u p (x) and energy huj p , and where n(p,t) 
j dx u p {x) k(x, t). 

We can then write the output energy flux in terms of 
the two-time correlation of the system. 



dt 



|^ = 2 \R{p)\ 2 »(jT*d« e--^-'V(i)a( U )>\ 



(11) 



This assumes that at time t = 0, the output field was in 
the vacuum state. 

When the output field is in free space and the only 
term in H Q is the kinetic energy, then the eigenstates 
are the momentum eigenstates. In this case, R(p, t) is 
just the Fourier transform of n(x,t). When there is a 
gravitational field, the eigenstates are Airy functions with 
a displacement which depends on the energy: 



u-(x) =Mk\\J3{x- ^)] 
mg 



(12) 



where Af is a normalisation constant, and the length scale 
is given by (3 — {2m 2 g /h 2 ) 1 ^ . In this case n{p,t) must 
be calculated numerically. 

We will now derive equations of motion for the system 
in the absence of pumping. 



III. SOLUTION IN THE ABSENCE OF PUMPING 

When the cavity is not connected to the pump reser- 
voir, we can find the Heisenberg equations of motion 
for the system operator, and then use Eq.(^) to replace 
il)H{x,t) with system operators: 



i(t) = —iu) a{t) + I dx k(x, t) x 



ipj(x,t) — / ds F(x,t, s)a(s) 



-iw a(t)- ds f(t,s)a(s))+£{t) (13) 

Jto 
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where 



f(t,s)— / dx n(x,t)F(x,t, s). 



(14) 



For most physical situations, the memory functions 
F(x, t, s) and f(t, s) are simply functions of (t — s), which 
means that Eq.@ is a Volterra equation of convolution 
type and can be solved using Laplace transforms if the 
Laplace transform of / exists. This has been solved ana- 
lytically for simple forms of the coupling in earlier work 
p8|^9| . In the limit where the decay of the system opera- 
tor a(t) becomes slow compared to the decay of the mem- 
ory function, the system operator can be taken out of the 
integral. This is a Markov approximation, and leads to 
exponential decay of the trapped mode. The general so- 
lution approaches the Markov solution as the strength 
of the coupling is turned down. Coupling the atoms out 
at a finite rate causes a significantly non-exponential de- 
cay. The difference between the exact and the Markov 
solutions is examined in detail in parallel work |po[ . 

The most surprising aspect of the exact solution is that 
it does not necessarily go to zero in the long time limit. 
If we continuously couple a single-mode in a trap to free 
space, then some of the initial population will disperse, 
but a certain fraction of the initial trapped atoms will 
remain in a non-dispersing eigenstate of the system and 
the coupling. Even if a momentum kick is given to the 
atoms as they are coupled out, the bound eigenstate ex- 
ists, although in this case it will have a lower occupation. 

The mechanism for the production of this bound eigen- 
state can be best seen in a dressed state picture. In this 
picture, the mixed state which was originally in free space 
experiences an attractive quasipotential. Let us find the 
eigenstate of this system explicitly: 

As all of the evolution is coherent and we are not in- 
cluding interatomic interactions, the atoms are acting in- 
dependently in first quantised picture, and we may use 
single particle quantum mechanics to describe the eigen- 
state. The state of the system is then described by a 
wavefunction 



I*) = ColVO + / dxij)(x)\x) 



(15) 



\lpE) = C E \lpa) + dx 1p E (x)\x) 



(17) 



The eigenvalue equation H\^je) 
the following equations: 



E\ipE) then leads to 



dy g*{y) 4>e{v) + c E huj Q = E c E , 



c E g(x) 



h 2 d 2 ^ E 
2M dx 2 



E ip E (x). 



(18a) 



(18b) 



These are most easily solved in Fourier space, so we 
take the Fourier transform of both sides and define the 
Fourier transforms: g(k) = 1/^/2/k j dx g(x) exp(— ikx), 

■0£;(fc) = l/v2~7r J dx ij)(x) exp(-ikx). We then show 
that the eigenvalue solves the equation 



E = TwOn 



2M 



dk 



\g(k)\ 2 

2ME 



k 2 - 



and the corresponding eigenfunction is 



2Mc ? 



m 



k 2 



2ME 



(19) 



(20) 



These equations can be solved analytically for some 
forms of the coupling g(k), and in general they can be 
solved numerically. For reasonable physical parameters, 
there is a single bound eigenstate with negative energy. 
As the coupling is reduced, this eigenstate becomes more 
weakly bound. If there was a gravitational potential in 
the output field, then we would expect that these bound 
states would become metastable, and would eventually 
decay. This can be verified numerically (2(]]. We would 
also expect atom-atom interactions to destroy this bound 
state, and in current experiments this is certainly the 
case, as the mean field interaction energy is much greater 
than the initial kinetic energy. 

In the following figure, we show that when grav- 
ity is added to the model, the bound state docs 
indeed decay. We have chosen a Gaussian cou- 
pling with width <7fc, which has the form g(x) = 
h^/T (2al/n) 1 / 4 exp(—(x (Jk) 2 )- We chose a weak grav- 
itational field so that the short time behaviour is un- 
changed. 



where \ip a ) is the lasing mode and ip(x) is the wavefunc- 
tion outside the trap in the position basis, {|x)}. 
The full Hamiltonian is 



h = nu! \ip a )(i(j a \ + — 



dx iff{x)\x)(i> a \+g*(x)\i> a )(x\) 



(16) 



where M is the mass of the atom and g(x) is the shape of 
the coupling. This Hamiltonian can describe a momen- 
tum kick on the atoms as they leave the trap by having 
a Fourier transform of g{x) which is not centred around 
zero. The trapped eigenstate will be of the form 
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FIG. 2. This figure shows the decay of the bound state 
when gravity is introduced. The solid line shows the evolu- 
tion of the number of atoms in the trap without gravity, and 
it approaches a steady state value. The dashed line shows 
how it decays in the presence of gravity (g = 9.8 sin(7r/20)). 
The trap atom number (a'(t)a(t)) is normalized so that 
( O t(0)o(0)> 



1. Parameters are F = 10 6 s 2 , = 10 6 m 1 , 



m = 5 x 10~*kg, L0 o = 2tt x 123 s" 



If we solve Eq.(18a) and Eq.(18b) for this example, we 



find that there is a bound eigenstate of energy E/h — 
-4.3 x lCPs" 1 with \c E \ 2 = 0.34. This means that 0.34 
of the original state of the system was in the eigenstate, 
and the rest has gone. It also means that 0.34 of that 
remaining eigenstate will be in the trap at any one time. 
This means that 0.34 2 = 0.12 of the initial population of 
the trap will remain in the trap in the long time limit, 
which is exactly what we see in our simulation. 

It is perhaps worth reiterating here that the analysis 
of this non-Markovian system is almost identical to that 
found in the field of photonic band gaps. Identical popu- 
lation trapping (in this case photons rather than atoms) 
has been calculated using identical analytical methods 
p2| , f23"|| . There are problems in using the Laplace trans- 
form method for slightly generalised models, however, 
even though the equations look almost identical. Particu- 
lar attention must be paid to the abscissa of convergence, 
as the Laplace transforms of many physically interesting 
kernels do not exist at all. A numerical method based 
on Laplace transforms will give incorrect answers in this 
case, and alternative methods must be employed. 

Gravity is not the only possibility for modelling an 
output coupling which will eventually remove all of the 
atoms from the trap. In practice, interatomic interac- 
tions may also destroy the stabilty of this state. The 
largest potential seen by the output atoms other than 
gravity will be a repulsion from the trapped atoms. This 
is because the output atoms will be relatively dilute com- 
pared to the atoms which remain in the trap. This can 
be modelled by an external potential for the atoms which 
will repel them from the trap. We will assume that both 
of these effects are small so that we can continue to use 
our single-mode approximation for the lasing state. 

The existence of a trapped state when there is no 
pumping will lead to concern later, as it will mean that 
pumping will cause the population of the trapped state to 
increase indefinitely, and no steady state will be reached. 

When the model is generalised to include the possibil- 
ity of pumping, the individual Heisenberg operators a(t) 
in the two-time correlation can no longer be calculated, 
so the same techniques for finding a solution do not apply. 
In the next section we will derive the equation describing 
the two-time correlation of the cavity field in the presence 
of pumping. 



IV. EQUATIONS OF MOTION INCLUDING 
PUMPING 

If the pump reservoir is sufficiently isolated from the 
cavity and external fields, and the pumping process is 
designed to be irreversible, then we may trace over the 
reservoir states to produce a master equation for a re- 
duced density matrix which describes only the cavity and 
the external fields. One example of a pumping process 
which can satisfy these requirements can be found in our 
earlier model, where cooled atoms in an excited state are 
passed over the trap containing the lasing mode JlJ] . The 
photon emission of the atoms would be stimulated by the 
presence of the highly occupied ground state, and they 
will make a transition into that state and emit a photon 
pO| . For a sufficiently optically thin sample, which can 
be made possible by having a very tight, effectively low 
dimensional trap, the photon is unlikely to be reabsorbed, 
and the process is effectively irreversible. We choose to 
model an optical cooling process rather than the more 
experimentally successful evaporative cooling process as 
we are particularly interested in designing a continuously 
pumped system with a steady state. 

After we have traced over the pump modes, we will pro- 
duce a term due to the effect of the pump in the master 
equation for the reduced density matrix p. It is impor- 
tant to note that we are not tracing over the output field 
modes, so our reduced density matrix spans the output 
field as well as the cavity field. 

We derive a pumping term based on an approach sim- 
ilar to that followed by Scully and Lamb |[l| and found 
in standard quantum optics texts [p9[ . We model pump- 
ing by the injection of a Poissonian sequence of excited 
atoms into the atom laser. These atoms may sponta- 
neously emit a photon and make a transition into the 
atom lasing mode. Alternatively, they may make a tran- 
sition into other modes of the lasing cavity. For simplic- 
ity, we consider an effective two-mode approximation. To 
obtain the pumping term, we consider the effect of a sin- 
gle atom injected into the atom laser, and then extend 
this to describe the effect of a distribution of atoms. This 
gives a master equation describing pumping in the num- 
ber state basis of the lasing mode p nm = (n\p\m), where 
each of these elements still include the output field. This 
equation can be expressed as 

(P)pump = (n. + Ala^y 1 p. (21) 

where the superoperators T> and A arc defined by 



V[c} = J[c]-A[c\, 
J[c]p = cpc\ 

A[c]p — -(c^cp + pc^c). 



(22) 
(23) 

(24) 



This is the same form as presented for a generic laser 
master equation by Wiseman p2[ where r is the rate 
at which atoms are injected into the cavity, and n s the 
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saturation boson number. In our particular model n s 
depends on the ratio of the probability that an atom 
will spontaneously emit into the lasing mode and the 
probability that the atom will emit into another mode. 
This pump process is Markovian for the same reasons 
that damping in optical lasers is Markovian, as there are 
photons being lost from the system. 

We may write this pumping term in the number basis, 
Pn,m = (n\p\m): 



(Pn,m) 



pump 



n s + (n + m) /2 
(n + ™ + 2)/2 
n s + (n + m + 2)/2 



Pn — 1 ,m—l 



' pn , m 



(25) 



At this stage, the elements p n ,m are operators which con- 
tain the state of the output modes. For optical lasers, we 
usually make the Born and Markov approximations at 
this stage. This is not possible for atom lasers except in 
extreme parameter limits, but we will review the results 
than can be obtained with these approximations so that 
we can compare them with the exact results. 



Pn-1 



r(n + 1) 
n + n s + 1 

ss 



(28) 



n-fBM P 



„ss 



+(n+l)7j££ 1) 



where 



lBM = 2 du R(/(u)e*""') 
Jo 



(29) 



which gives us the steady state for the atom number dis- 
tribution in the cavity: 



n SS_ Af ( r hBM)' 



(n + n s )\ 



(30) 



where Af is a normalisation constant. This form is al- 
most identical to that obtained for the optical laser |29] . 
The distribution looks thermal for a = r/7 < n s , and 
for a > n s , which is called the above threshold regime, 
it approaches a Poissonian distribution with mean atom 
number n and variance V given by: 



A. Pumping using the Born approximation 

The Born approximation assumes that the density ma- 
trix p can be separated into the product of the reduced 
density matrix for the system a, and the reduced den- 
sity matrix for the output field. It further assumes that 
the output field remains in its original state. If we then 
trace over the output field modes, we obtain a master 
equation for the lasing mode which has the same pump- 
ing term as given above p0| . The reduced density now 
describes only the cavity mode, and a n , m — (n\a\m) are 
c-numbers. The damping term has the form: 

Wdamp = ~ ! du [f(u)e^(a^aa(t - u) 



-aait — u)a') + h.c. 



(26) 



where f(t) has been defined in Eq.(|Tj). 

If we combine the pumping and damping terms, we can 
generate equations of motion for the steady-state atom 
intracavity number, p n = <r nn : 



rn / x r ( n + 1) f \ 

P „_!(t) rrPn(t) (27) 



n + n 



n + n s + 1 



-2n / dw5i(/(u)e M °>„(i-u) 



Pn(t) 



+2(n + l) / d M 5R(/(u)e M ° £ K +1 (<- U ). 
Jo 

where 3i(g) is the real part g. 

By setting the derivatives to zero and assuming that 
the functions p n {i) approach a constant p^ s as t goes to 
infinity, we obtain a recursion relation 



IBM 

V = n + n s 



(31) 



Eq. (^) will give us the properties of the output field if 
we can calculate the two-time correlation of the intracav- 
ity field. We proceed by finding the equation of motion 
for the expectation value of the field operator: 



dt 



pump n=0 

00 

= ^ Vnpn-l,n 



r 



2(n + n s ) + 1 



2n + 2n s + 1 

(32) 



where we have used the fact that the number distribution 
is localised to replace n by n in the denominator. By 
expanding the exact expression as a Taylor series in 
we see that there is an error term of the order of 1/n 2 . 
This is awkward, as if the first order terms cancel then 
the lincwidth will be of this order. This is in fact what we 
are hoping to discover, as it will mean that the spectrum 
will become more narrow as we increase the pumping of 
the laser. If this occurs, then we may approximate the 
magnitude of the off-diagonal elements of p as though it 
were a coherent state. We then find that the expression 
given in Eq. (^||) is correct to third order. 

The analysis which produced this pumping term was 
only made possible by assuming that the cavity atom 
number becomes localised around some mean. If we do 
not make the Born approximation, then this fact must 
be taken on faith, and later shown to be a self-consistent 



G 



solution. Furthermore, if we do not make the Born ap- 
proximation, then we do not know at this stage what the 
mean value n will be in terms of the other parameters. 
Once we accept that the atom number distribution peaks 
around a particular value, however, then we may calcu- 
late the output spectrum without having to trace over 
the output modes. 



B. Linewidth using the Markov approximation 

The Markov approximation can be made in certain 
physical limits |20| . It corresponds to assuming that the 
decay of the reservoir correlation function f(t) is fast 
compared to the rate of change of the density matrix. We 
wish to describe the atom laser in the parameter regimes 
where it becomes invalid but we will compare the results 
of this simple calculation to the results of a calculation 
which does not make this approxiamtion. 

If we make the Born approximation, trace over the out- 
put modes and then further assume that a Markov ap- 
proximation can be made, then we can write the damping 
term of the master equation as: 



(*)d 



amp 



7bm£> [a] a. 



(33) 



The derivation of this equation has been examined in de- 
tail in other work PQl , and it has been assumed to be of 
this form in earlier atom laser models. It can be gener- 
ated from Eq. (p6|) by assuming that u(t — u) = cr(t), and 
can therefore be taken out of the integral. The damping 
constant Jem is the same as the one defined in Eq.@. 

The total equation of motion for the expectation value 
becomes 



dt 



2(n- 
r 

An 2 



,t\ 



2fl) (at) 



(34) 



which has an error term proportional to r/n 3 . 

The solution to this equation is an exponential decay, 
and the energy spectrum can therefore be shown to be 
Lorentzian, with a width of T = r/(4n ), which is of 
the order of 7bm/^j and will therefore go to zero as we 
increase the pumping of the laser. This is called gain- 
narrowing, and it is a well-recognised feature of the op- 
tical laser, but even in optical lasers it is a feature which 
only exists for certain pumping models. If the pump 
has a significant response time, then an optical laser can 
have a linewidth which scales as r/n, which is simply 
proportional to the cavity linewidth ^9|. Alternatively, 
a well-designed pump can actually lead to subPoissonian 
statistics |29[| . 

Since this model of the atom laser exhibits gain- 
narrowing when we make the Born-Markov approxima- 
tion, we hope to find that it also exhibits gain-narrowing 
when we solve the model correctly, as it is this feature 



which allows the atom laser to p roduce output with a 
high spectral density. In Sec. VII we shall show that this 
is indeed the case. 



C. Correct treatment of pumping 

If we do not make the Born or Markov approximations, 
but we do assume that the trap population is localised 
around some (at this stage unknown) value n, then we 
can use Eq. (p2[) to produce a more general equation of 



motion for the expectation value 



»t\. 



dt 



(ot)(t) = (ta; + P)(ot)(t)+ Ida K*(x,t)(^){x,t) (35) 



where 



P 



2(n + n s ) + 1 



(36) 



well above threshold. Remember that we can no longer 
relate n directly to the physical parameters of the prob- 
lem using Eq.(|3l|), which has used the Born approxima- 
tion. If the solution exhibits gain-narrowing, then the 
value of P that we use will actually determine our value of 
n s . We are therefore producing a solution method which 
can be solved iteratively to produce a self-consistent so- 
lution rather than directly generating it. 

Under our assumptions the pumping is effectively lin- 
ear, so we may use the quantum regression theorem. We 
then recall Eq.(J^), and we may immediately derive the 
following integro-differential equation for the two-time 
correlation: 



£<at(t + r)a(t)) = {ilu + P)(ct(t + r)a(t)) 

OT 

- J Q t+T duf*(t + T- u)( a t(«)o(t)) 



(37) 



where r > 0, and f(t) has been defined in Eq.(jl4|). We 
have set to to zero, and assumed that at this point in time 
there were no atoms in the output field. In other words, 
at t = the coupling was switched on, and the expec- 
tation values involving the normally ordered operators 
tpi(x,t) and tpj(x,t) will be zero. 

This equation is not sufficient to specify the dynam- 
ics of the cavity, as it is only a single partial integro- 
differential equation in a two dimensional space. We also 
require the integro-differential equation for the intracav- 
ity number, which we can generate in a similar manner. 
Well above threshold, we obtain: 



d_ 

dt 



; a to)(<) 



du2Re{J*{t~u)(a\u)a(t))} (38) 



These equations are difficult to solve in general, but 
can be solved simply in various limits. For example, if 
the kernel f(t) was a 5-function then the equations would 
become local and the solution would be a simple exponen- 
tial. This is the case for the Markovian example shown 
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in Sec. IV B above. In the broadband limit of the opti- 



cal theory, the function f(t) is equivalent to the Fourier 
transform of a constant, which is exactly a <5-function. 
Although the broadband limit can be a good approxima- 
tion for the atomic case as well [Q, the fact that atoms 
will disperse in free space means that the system has an 
irreducible memory, and is non-Markovian. 



D. Memory functions 

Our model for the output involves a change of elec- 
tronic state via a Raman transition, which has recently 
been achieved experimentally ]l7| . From a specific model 
of this coupling we can derive the memory functions 
f(t,s) and F(x,t,s). The form of f(t,s) in this case 
has also been derived and discussed in the recent work 
by Jack et al. p3fl . In position space, the shape of the 
coupling will be defined by the shape of the lasers mul- 
tiplied by an envelope of the spatial wavefunction of the 
trapped state. For simplicity of calculation, we will as- 
sume that the coupling is Gaussian in form, and that 
there is no net momentum kick given to the atoms. This 
will allow us to produce analytical forms for the memory 
functions. Let the coupling be defined by 



_/2a 2 \ 1/4 
K i x t t) = \/7 [ I < 



(39) 



where Turk is the momentum width of the coupling and 
7 is the strength of the coupling. 

In the presence of a gravitational field, V = mgx, the 
Green's function in Eq.(|8|) can be found as a standard 
result §3: 



G(x,t,y,s) 



1 



Aiti\(t - s) 

f i(x-y) 2 ig(t- 3 )(x+«) ig- ■■ 
p} 4A(t-s) 4A 



(40) 



2 (t-s ^ 3 - 
ISA 



where A = Ti/(2M). 

This leads to the following forms for the memory func- 
tions: 



F{x,t,s) = 

and 
/(*,*) = 



o\V 4 

2 \ / 7<T fc 



x 



iv j y 1 - 4 i ak A At 

At g(At 2 g + 6 x)+i (At 4 g 2 + 12 At 2 g x-12 x 2 ) A 

e 



(41) 



12 A (i+4 At A err' ) 



7 exp(-^) op(^L 

k 

y/l + 2 i A At al 



(42) 



where At = t — s. 

If we make g = 0, then we have the memory functions 
for coupling into free space. We can see that /(t) will 
then go as 1/y/t in the long time limit. The broadband 



limit may also be found by taking <Jk — > oo while j/fTk 
goes to a constant. This limiting case has been examined 
in detail in previous work [ ^8|JlS| ] . For a non-pumped cav- 
ity with a broadband output coupling, the long time limit 
of the output spectrum has three regimes. In the limit 
where the coupling is much faster than the free field dis- 
persion rate of the atoms, the energy spectrum looks like 
the cavity wavefunction which has appeared in free space. 
In the limit of coupling that is very slow compared to the 
dispersion rate, the output approaches a Lorentzian. As 
the coupling strength is increased from one limit to the 
other, there is a reasonably complex behaviour which at 
first looks like a deformed Lorentzian, and then produces 
a fringe-like structure in the spectrum. Eventually, the 
oscillations decay and the spectrum begins to look like 
the cavity wavefunction. 

There are no useful approximations which can simplify 
the full calculation. For the purposes of numerical calcu- 
lations when g^O, using the broadband limit is actually 
less tractable than using the more realistic memory func- 
tion. This is because the integrals involving F(x, t, s) and 
/(t, s) become unbounded in amplitude, and their con- 
vergence is due to their highly oscillatory nature. Retain- 
ing the full form of the equations includes the Gaussian 
envelope, which defines a natural upper bound to the in- 
tegrals. Since the oscillations grow rapid on the same 
scale as this envelope decays, it is probable that the cor- 
rect solution and the broadband limit are qualitatively 
very similar. 

If the two-time correlation is very slowly changing then 
it may be removed from the integral in Eq. ( |37j ) , which is 
equivalent to approximating the kernel to be extremely 
narrow. We may therefore guess that well above thresh- 
old, where we hope to find a slowly decaying two-time 
correlation (or a narrow linewidth), the Markov approx- 
imation may hold, and the output will be Lorentzian in 
energy. However, there is also the possibility of other 
behaviour. Even if the amplitude of the two-time corre- 
lation decayed slowly and exponentially, a fast change in 
the phase of the solution would not only mean that there 
was a shift in the resonance, but would also significantly 
change the effective decay constant. There is no obvious 
way to determine whether the coupling will induce such 
a rotation without actually solving the original Volterra 
equation. In fact, we must solve Eq. ( |37| ) and Eq. (S3) 
self-consistently to be sure that this system will exhibit 
gain-narrowing of the output field. 

In the next section, we will show how we may find a 
solution to Eq. (37) and Eq. (j38|), and then in the fol- 
lowing sections we will calculate the output from a laser 
with g — 0, and then with a nonzero gravitational field. 



V. SOLUTION METHODS 

Equation (p7|) and Eq. (|38| ) do not form a standard 
pair of partial integro-diffcrcntial equations. The deriva- 
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tive in Eq. ( |37| ) is only defined for r > 0, and so we cannot 
require that the solution obey this equation through the 
whole domain of the integral. This means that Eqn. j37| ) 
cannot be integrated to find the solution, as we do not 
actually know the derivative at any point. Since the two- 
time correlation is Hermitian, we can we rewrite the in- 
tegral so that the domain remains in the r > plane, 
but we still do not have a continuously defined derivative 
along the length of the integral. This is shown graphi- 
cally in Fig. (0). 




FIG. 3. The r > half plane of the two-time correlation. 
Equation (j37j) defines the horizontal partial derivatives over 
the entire region, and the domain of the integral correspond- 
ing to the derivative at point P is shown with a dashed line. 
The only derivative known at point P' is in the direction 
shown. Equation (^) defines the derivative along the diag- 
onal, and the domain of the integral corresponding to the 
derivative at point Q is shown with a dotted line. 

We proceed by making an ansatz which uses the solu- 
tion of Eq. (|37]) which has been extended into the region 
r < 0. We then use the r > portion of this solution to 
substitute into the two-time correlation in Eq. (^8|). We 
introduce the function J(t), which is the solution of the 
equation 

^ = (iuj + P)J(t) - J dw f*(t, w)J(w). (43) 
This means that 



(a f (t + i>(f)) = (a f a)(t)^ 



Jit) 



(44) 



is a solution of Eq. (37) with the correct initial condition 
at t = 0. We then substitute this result into Eq.(^8|): 

<it =r ~ Jo rfw(aa>w H V) > (45) 

Solving these two equations gives the two-time corre- 
lation for the lasing mode, from which we may find the 
properties of the output field. It is only consistent with 
our linearisation of the pumping if the number of atoms 
in the trap, n, is stable around the value which originally 
produced the parameter P. Since we require P to gen- 
erate the solution, and n is simply the long time limit of 
Eq. (f45|), the effective free parameter is n s . This thresh- 
old parameter must be much smaller than n, so we search 
for a value of P which gives the result h^> n s . 



Once it is established that Eq. (45) is approaching a 
stable steady state, a fast way of finding that steady 
state is to set the derivative to zero, and assume that 
(a*a)(w) — h over the support of the kernel. This gives 

(46) 



There are many numerical methods for sol ving Volterra 
equations, and since f(t, s) — f(t — s), Eq. ( |37|) is of the 
convolution type, which has many specialised methods of 
solution. We have found that this equation can be solved 
analytically by a Laplace transform method in the ab- 
sence of gravity and in the broadband limit. However, 
for the more general forms of coupling that we have ex- 
amined in this paper, even numerical Laplace transform 
methods become invalid as the Laplace transform of the 
memory function does not always exist. In these cases 
we require an alternative numerical method. 

We transform Eq. ( f43| ) into a rotating frame by intro- 
ducing the function 



N(t) = J(t) exp(-iw t - Pi) 



(47) 



and rewriting the equation of motion for J(t) in terms of 
N(t): 



dN(t) 
dt 



dw H(t,w)N(w) 



(48) 



where H(t,w) = f*(t, w) exp(— [iu + P](t — w)). Since 
f{t,w) is actually just a function of (t — w), this equa- 
tion is of convolution type. Unfortunately, the numeri- 
cal methods for finding the solution of convolution type 
Volterra equations depend on the simplicity of the kernel 
for their effectiveness. We have found that the most effec- 
tive numerical method for solving this equation is simply 
a direct integration of Eq. (48) using a second order al- 
gorithm for both the integration and the calculation of 
the integral to find the derivative at each timestep. Such 
techniques can be found in many common collections of 
numerical methods (jj5|. 

In the following section, we will examine the case where 
g = 0, and atoms are simply diffusing away from the trap. 



VI. IDEAL ATOM LASERS MUST POINT DOWN 

Real atom lasers will be able to point in any direction 
and still expect a beam of atoms to be emitted. The 
three reasons for this are repulsive atom-atom interac- 
tions which will eject the unconfined atoms, the presence 
of gravity which will accelerate the atoms away from the 
lasing mode, and any momentum kick given to the atoms 
by the coupling process. We shall show that it is not actu- 
ally possible to produce even an idealised model without 
either the gravity or the atom-atom interactions. 

We will do this by attempting to make a model of an 
ideal atom laser without considering complicating factors 
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such as atom-atom interactions and gravity. In this limit, 
we can use analytical methods to solve Eq (^7|) . 

If we take the Laplace transform of Eq. (|37D, and then 
use the convolution and derivative theorems, then we ob- 
tain 



j(t) = j(o) cr x 



1 



s -(iu + P) + £{/*} (s) 



(*) (49) 



where £{/}(s) = f^°dt /(i)exp(— st) is the Laplace 
transform of f(t), and C~ l is the inverse Laplace trans- 
form. 

When .9 = and we go to the broadband limit, then 
the memory function is 



i 



(50) 



where fHs a constant related to the strength of the cou- 
pling. 

For this form of f(t), we can find the solution to 
Eq. (^) using results from eariler work |lj||2(|. We dis- 
cover that in the long time limit it grows exponentially. 
The number of atoms in the cavity also grows exponen- 
tially when this solution is substituted into Eq. ( |38| ) . This 
shows that for g = 0, the steady state approximation 
on which we based our linearisation of the pumping in 
Sec. (IV C) is not self- consistent. 

The exponential growth is not physically realistic, and 
is merely an artifact of the breakdown of our approxi- 
mation. In essence, our approximation will ignore the 
depletion of the pump if the linearisation is not valid and 
so in that case we would be modelling a pump which can 
deliver an infinite atomic flux. The reason for the ab- 
sence of the steady state is clear once we have discovered 
the existence of the trapped state which was described in 



Sec. (III). Since our trap does not empty when there no 



pumping, it is not unreasonable to expect that a certain 
fraction of the atoms coming into the trap from the pump 
will enter the trapped state. This will mean that the 
number of atoms in the cavity would continually grow. 
We expect that solving this system without linearising 
the pump would lead to a solution that involved a trap 
number which increased linearly. 

The absence of a steady state means that this model is 
not suitable for describing an atom laser, and a greater 
level of detail is required in order for the model to be- 
have realistically. In other words, our idealised atom laser 
model must be more complete than the idealised optical 
laser model. We also realise that it is the nature of the 
output coupling which must be made more realistic. A 
purely damped cavity must lose all of its atoms in the 
long-time limit, or there is not really any hope that the 
pumped atom laser will reach a steady state. 



Our examination of the trapped state in Sec. (HI) 
guides us here, as we know that merely adding a momen- 
tum kick to the coupling does not remove the trapped 
state. It is likely that appropriately strong and repulsive 



interatomic interactions would provide the required real- 
ism, but we find that that the effects of a gravitational 
field are simpler and easier to model. In the next section, 
we show that an atom laser in a gravitational field does 
indeed produce a self-consistent steady state. 



VII. OUTPUT FROM AN ATOM LASER 

The output from an atom laser where the output atoms 
accelerate under a uniform gravitational field is a func- 
tion of many physical parameters. In terms of indepen- 
dent timescales for the dynamics, we can control the dis- 
persion rate in free space (by choosing the mass of the 
atoms); the time taken for the output atoms to leave the 
trap due to gravity; the coupling rate between the trap 
and the output; the pumping rate; the threshold pump- 
ing rate and the trap frequency. Within this phase space, 
we are constrained by physical limits such as the choice 
of atom and atomic trap. We are also constrained by our 
operational limits, such as our maximum allowable trap 
density and our desire to be well above threshold. Fi- 
nally, we are constrained by computational limits, which 
make it difficult to calculate the behaviour of the laser 
when any of these timescales become radically different 
from the others. In principle it is possible to discover 
approximations that will work in any of these limits, but 
each one would have to be considered separately. We will 
simply examine the behaviour of the output for many dif- 
ferent pumping and damping rates. We will choose the 
remaining parameters as realistically as possible. 

The shape of the coupling n(x,t) is determined largely 
by the spatial wavefunction of the laser mode, and this 
in turn is a function of the trap energy. We use the 
trap frequency lo = 2tt x 123 Hz, which is near those 
used in experiments . We use an atomic mass of to = 
5x 10~ 26 kg and a gravitational field of g = 9.8 sin(0.18) m 
s~ 2 . We assume the coupling to be Gaussian, with a form 
given by Eq. (|39|), and a momentum width = 4.4 x 10 5 
m _1 . 

Our first calculation will show that the system can 
reach a steady state. We choose a pumping rate of 
r = 5.0 x 10 s -1 , and a threshold cavity number of 
n s = 8. The damping rate is then chosen so that the 
number of atoms in the trap will be much larger than n s . 
We use 7 = 2.0xl0 5 s~ 2 , which gives n = 1180 > n s . We 
then solve Eq. ( p8| ) and use the result to solve Eq. (|3^), 
which gives us the number of atoms in the trap as a func- 
tion of time. The results are shown in Fig. (^), where we 
can see that after initial transient behaviour, the number 
of atoms reaches a steady state. 
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FIG. 4. This figure shows that the pumped atom laser 
reaches a steady state when the output field experiences a 
gravitational field. The solid and dashed lines indicate the 
solution for different initial trap numbers. 

The initial transient behaviour is only correct as a first 
approximation, as we are actually determining a self- 
consistent steady state for this system in which the pa- 
rameter P is stable. The different starting values for n 
are included to demonstrate that this solution can be 
found in a stable fashion. 

The two-time correlation derived in this calculation is 
approaching an exponential in the long time limit. Unlike 
the Born-Markov solution, this exponential has a large 
rotating term. This rotation gives the output a frequency 
shift and it also alters the integral in Eq. (^) which in 
turn affects the rate of decay. Even in the strong pumping 
limit, where the amplitude of the two-time correlation is 
decaying very slowly, the Markov approximation gives an 
incorrect result. 

In a counterintuitive way, the fact that the Markov 
approximation is not valid makes it harder to compare 
the correct calculation to a Born-Markov result. When 
the memory function is effectively a delta function, then 
the Markov approximation gives the same result in any 
rotating frame, but when the memory function doesn't 
decay instantaneously, then the decay rate and the fre- 
quency shift depend on the rotating frame in which the 
Markov approximation is made. This is nothing more 
than saying that the Markov approximation is not self- 
consistent. For the purpose of being definite, we will use 
the most "natural" version of a Markov approximation, 
in which we rotate at the syst em fr equency to . This gives 
the linewidth shown in Sec. [VB and the damping rate 
given by the Born approximation, Eq. (^9|). 

We now look at what happens when the pumping is 
increased. We use a damping constant of 7 = 2.0 x 10 4 
s -2 , and a threshold of n s — 47, and vary the pumping 
rate. As the pumping rate increases, the steady state 
number of atoms in the trap increases, and the modu- 
lus of the two-time correlation decays more slowly. The 



energy spectrum of the output flux is proportional to 
the Fourier transform of the two-time correlation through 
Eq. (|ll]), so we can see that the energy spectrum is be- 
coming more narrow as the two-time correlation becomes 
broader. This means that the laser is experiencing gain- 
narrowing. 

In Table 1 we show the results of these calculations. We 
give the pumping rates r and the resulting mean atom 
cavity numbers n. The linewidth of the output energy 
flux, r, is calculated directly from the two-time correla- 
tion using Eq. (pd|). This is compared to the linewidth 
given by the Born-Markov approximation Tbm, which 
was found from Eq. ([34|). The spectral shift Au; is the 
amount by which the correct result is shifted from the 
Born-Markov result. This spectral shift is seen to be 
largely a function of the damping only, and therefore does 
not change as we vary the pumping. 

We claimed in the introduction that a continuously 
pumped atom laser can have a linewidth much narrower 
than can be obtained by dropping the atoms from the 
trap after a rapid state change. In the fast coupling limit 
the output field will simply look like the original con- 
densed wavefunction. More precisely, the output field 
will have the same shape as n{x), which for our example 
will be a Gaussian, given by Eq. (|39|). We may proceed 
to find the energy spectrum by changing to the free space 
energy basis (uj — \k 2 ). The resulting energy spectrum 
is 



1 



ex P(^T32)- 



(51) 



which is normalised so that J °° duo \^{w)\^ = 1. 

This spectrum is singular, so a FWHM definition of 
the width would be meaningless. We shall define the 
linewidth, r^- ast , obtained by a rapid state change by the 

equation J Q r/ ° st duj \if;(w)\ 2 = 1/2. For the parameters 
used above, this gives us r/ ast = 88 s . Comparing this 
linewidth to those calculated in Table 1 shows that our 
atom laser gives an improvement of two to three orders 
of magnitude in spectral density. 

TABLE I. Linewidth as a function of pumping 



r (10 3 /s) 


n 


r(s- A ) 


r SM (s- 1 ) 


Alo (s" 1 ) 


20 


450 


2.1 


0.025 


0.16 


40 


910 


1.1 


0.012 


0.15 


80 


1800 


0.56 


0.0062 


0.14 


800 


1.8 x 10 4 


0.035 


0.00062 


0.14 
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We plot the spectral flux corresponding to three of 
these pumping rates in Fig. The vertical scale is 

normalised to the peak height for each plot so that the 
width of the spectra can be easily compared. 
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FIG. 6. The output energy flux spectrum without (solid 
line) and with (dashed line) the Born-Markov approximation. 
Parameters are r = 5.0x 10 5 s~ 1 , 7 = 2.0x 10 5 s~ 2 and n s = 8. 



FIG. 5. The spectra of the output energy flux for three 
different pumping rates. The solid, dashed and dotted lines 
represent r = 800 s" 1 , r — 80 s -1 and r = 40 s _1 respectively. 
Other parameters are given in Table 1. 

Although this line narrowing seems to suggest that the 
linewidth of the output can be reduced indefinitely, in 
practice there will be a limit due to technical noise in 
the pumping, output coupling and the laser mode. A 
limit due to the finite temperature of the trap has been 
calculated recently by Graham j37j. 

From Table 1, we can see that the solution obtained 
with the Born-Markov approximation is significantly dif- 
ferent from the correct solution. In Fig. (||) we plot 
the output energy spectrum with and without the Born- 
Markov approximation. We have chosen the parameters 
r = 5.0 x 10 5 s _1 , 7 = 2.0 x 10 5 s" 2 and n s = 8. This 
example is well above threshold, with a mean atom num- 
ber of n = 1200. We scale the two solutions so that 
their peaks are the same height so that they can easily 
be compared. 



VIII. CONCLUSIONS 

We have demonstrated that for some interesting pa- 
rameter regimes an atom laser cannot be modelled by 
Markovian equations. Although this was shown ear- 
lier for a nonpumped cavity |lq-|20|, there was some 
hope that the slowly decaying two-time correlation in the 
pumped case would resurrect some form of the Markov 
approximation. We have shown that this can only be 
true when the coupling is very weak, which will tend to 
make either the atomic flux very weak, or else require 
extremely large atomic densities. 

For noninteracting atoms, the non-Markovian equa- 
tions can be solved analytically in the limit of zero pump- 
ing, as was the case for the pulsed atom laser p5[ . We 
have developed a numerical method for calculating the 
output spectrum in the presence of pumping. We have 
found that when the full, non-Markovian dynamics are 
considered, it is necessary to add gravity to the model in 
order for a self-consistent steady state to be produced. 

This model does not include atom-atom interactions 
and therefore only works when the atomic field is very 
dilute, which is consistent with our parameters. The 
advantage of this model is that we have calculated the 
full quantum statistics of the lasing mode, which is im- 
possible with atom laser models which involve nonlinear 
Schrodinger equations based on mean-field theory. Fu- 
ture work will involve generalising this atom laser model 
to include some of the effects of interactions without mak- 
ing an initial mean-field approximation. This will allow 
us to determine the limitations on the coherence of a 
practical atom laser, and thus the limits on their inter- 
ferometric applications. 
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APPENDIX A: DERIVING THE LANGEVIN 
EQUATION 

In this appendix, we derive Eq. (Q). This means that 
we need to find tpu(x,t), which can be written: 

il> H (x,t) = U} nt {t,t ) V/OMo) U int (t,t ). (Al) 

where Ui n t(t,to) is the unitary evolution operator cor- 
responding to the interaction Hamiltonian, Eq (Q). It is 
the identity operator when t = to, and obeys the dynamic 
equation 



dU int (t, t ) 
dt 



Hi n t(t) Ui n t(t, to)- 



(A2) 



We proceed in a very similar manner to a calculation 
made by Jack et al. in their work on non-Markovian 
quantum trajectories ]3q]. Let us consider the operator 



U(t,t ) = T a {V + V V.}(t > t ) 



(A3) 



where T a is a time ordering operator which denotes a 
time ordering on the a operators only, and where we have 
defined the operators 



V + (t,t )=exp{- Jf £t( s ) aj ( s )j (A4) 

flt du /to dv f{u ' v) '' WaiM (A5) 
ds£(s)a\(s)) (A6) 



V 

V- = exp 



t 



where f(u, v) = [£(zt), £ '(v)] is the same function defined 
in Eq. @. 

Now {/(to, to) is clearly the identity operator. The 
equation of motion for U(t,to) is given by 



T a I V+Vo 



dV- 

-T a {tf(t)a r (t)V + V V-} 

-T a jv+ J ds f(t, s) 4(t) ai (s) V Q V- 

+T a {-V + V ^(t) ai (t)V-} 



T a {[V+(t,t ),£(t)]} = / duf(t,u)ai(u)V + (t,t ). 

Jt 

(A8) 

We therefore see that since U(t,to) obeys the same 
equation of motion and has the same initial condition 
as Ui n t(t,to), then it must be the same operator. We 
use this alternate form of the evolution operator to find 
ipH(x, t). 

t/> B (x,t) = UHt,t Q ) Tpi(x,t) U(t,t ) 

= ^(x, t) + U\t, to) T a {[^(x, t),V+] V V-} 
= 4> I (x 1 t) + uHt,t ) x 



dsF(x ) t,s)T a {a I (s)V + V V^} (A9) 



where F(x, t, s) = [ipi(x, t),£'(s)] is the same function as 
defined in Eq (j^). We then use the lemma 

T a {a z (s) V+VoV-} =T a {a I (s)V+(t,s)V (t,s)V-(t,s) x 

V + (s,to)V (s,to)V-(s,t )} 
= U(t,s) ai(s)U(s,t ) (A10) 



to reach the final form of Eq.(A9) 



ijiH(x,t) = ipi(x,t) — / ds F(x,t,s) Off(s). (All) 

J tg 



■j- H int (t)T a {V+V V-}(t,t ) 



(A7) 
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